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Abstract 

MnsGesCa; films with x > 0.5 were experimentally shown to exhibit a strongly enhanced Curie 
temperature Tq compared to Mn^Ge^. In this letter we present the results of our first principles 
calculations within Green's function approach, focusing on the effect of carbon doping on the 
electronic and magnetic properties of the MnsGes. The calculated exchange coupling constants 
revealed an enhancement of the ferromagnetic Mn-Mn interactions mediated by carbon. The 
essentially increased Tq in MnsGesC is well reproduced in our Monte Carlo simulations and together 
with the decrease of the total magnetisation is found to be predominantly of an electronic nature. 
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Intense efforts have recently been devoted to the fabrication of high-Tc metal/metalloid 
materials for their potential incorporation into spintronic devices. Among them, MnsGea 
seems to be a promising candidate due to its compatibility with mainstream silicon tech- 
nologyi. The low Curie temperature {Tq = 304 K-) of this compound, which is a main 
disadvantage for technological applications, has recently been overcome by carbon doping in 
MnsGeaCj; films.- It was shown^ that carbon incorporation into the octahedral voids of the 
hexagonal MusGes cell leads to a continous increase of Tc with the doping level in MusGesCa; 
and suggests a saturation when all octahedral voids are filled by carbon. A maximum Tc 
of 445 K was obtained for MnsGeaCo.s films, obtained by magnetron sputtering. The car- 
bon implanted MnsGeaCo.s films were found to exhibit magnetic properties very similar to 
their sputtered counterparts^. The average saturated moment of 2.2 /ib/Mu was observed 
for C-implanted MusGesCo.s films (1.1 /ib/Mu in the sputtered samples) and turned to be 
somewhat smaller than 2.6 /xb/Mu of MusGes polycrystals^. At the same time, the isostruc- 
tural antiferromagnetic MnsSia with Tn=98 K was reported to exhibit ferromagnetism when 
doped with carbon, with the transition temperature reaching 350 K in MnsSiaCo.s^. 

While in MnsSisCa: films the interstitial carbon leads to a lattice expansion, the increase 
of Tc upon doping in MusGesC^: is accompanied by a lattice compression^. However, it 
was suggested that the observed lattice distortions in both MusGesCa; and MnsSisC^: can 
have a very small influence on the increase of the transition temperatures^. On the other 
hand, it was also concluded that the double-exchange mechanism plays only a minor role in 
enhancing the ferromagnetism in MnsGeaCj;. In spite of the existing experimental data on 
the magnetic properties of C-doped Mn-Ge and Mn-Si compounds, the information about 
the local magnetic moments, magnetic interactions and the role of carbon interstitials on 
the enhancement of the Tq of MnsGesC^; is still lacking. Hence, in this letter we explore 
possible zero- and finite-temperature magnetic properties and transition temperatures of 
the pure and C-doped MnsGea compounds by using first-principles calculations and Monte 
Carlo simulations. 

We follow a standard scheme for the calculation of thermodynamical properties of mag- 
netic systems: (1) first the Green's function formalism and the magnetic force theorent^ are 
employed to determine the exchange integrals Jij from first principles; (2) these are used as 
input to a classical Heisenberg Hamiltonian of the form H = — J2i,j-i^j Jij^iej, where i and 
j are the site indexes and is a unit vector along a spin moment at i. The Tc is estimated 
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from the peak in susceptibility-temperature curve, obtained by Monte Carlo simulations 
using 2160 magnetic- atom supercells. The ab initio calculations are performed by the full- 
potential screened Korringa-Kohn-Rostoker (KKR) Green function method^ within the local 
spin density approximation^ of density functional theory. For structural relaxation we use 
the projector-augmented wave approach as implemented in the ABINIT code,- within the 
generalized gradient approximation for the exchange-correlation energy.— 

In order to study the effect of structure relaxations on the magnetic properties of MngGea 
and MnsGeaCj;, we performed calculations for both rigid and relaxed lattices. We have 
used the following lattice parameters for MusGes, as well as non-relaxed C-doped MnsGes: 
a = 7.184 A, c/a = 0.703^. The structure of the MnsGeaC^; system at x = 1.0 was relaxed at 
the experimentally obtained parameters of MnsGeaCo.rs films^: a = 7.135 A, c/a = 0.700. 
The unit cell of the hexagonal MusGes (space group PQ^/mcm, DSs-type structure — ) 
contains two sublattices of Mn (Mni and Mnn) with different coordination. The crystalline 
structure of the MusGcsC^. (Fig. [1]) was found to be similar to the MusGes structure with 
carbon atoms occupying interstitial positions at the center of Mnn octahedron^. 

According to neutron scattering experiments^, the difference in local environment of Mnj 
and Mnn atoms in MusGca is believed to be responsible for the different magnetic moments 
on them. The calculated spin moments on Mnj, Mnn and Ge atoms in relaxed (rigid) 
MnsGes are equal to 2.11(2.09), 3.11(3.15) and — 0.14(— 0.15) /ib correspondingly and are 
almost uneffected by structural relaxation. The smaller magnetic moment on Mnj is due 
to a direct Mnj-Mni interaction at a rather short distance (2.526 A).— The obtained values 
agree well with the previous calculationsiii^ and experimental observations^. The electronic 
density of states (DOS) (Fig. [2]) is metallic, in agreement with the recent calculations.- The 
partial DOS of Mn atoms is dominated by the 3(i-states, while the largest contribution to 
Ge DOS is due to 4p- and 4s- states. The DOS at the Fermi level (Ep) is dominated by 
the Mn d-states, pointing at the significant Mn-Mn interaction.- The exchange splitting is 
larger for Mnn ci-states compared to Mnj, in agreement with the larger magnetic moment 
on Mnn. 

In order to understand the role of C doping in the enhancement of Tq, we first calculated 
the structurally non-relaxed MusGesC unit cell with two carbon atoms in the octahedral 
voids and compared it to the corresponding results for the undoped compound. Such situa- 
tion represents the pure hybridization effect between the carbon and the neighboring Mnn 
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atoms. As it follows from Fig. [2l carbon incorporation into MnsGea essentially changes the 
3(i-states of Mnn, leaving the corresponding Mnj states almost unaffected. The hybridiza- 
tion between the C 2p- and Mnn Srf-states leads to a shift of the main occupied peaks in 
the majority, as well as unoccupied peaks in the minority 3(i DOS of Mnn towards E^-, 
consequently increasing the density of states N^E-p) at the E-p in both spin channels. The 
exchange splitting is now lower than in the pure MusGes, resulting in a reduced magnetic 
moment on Mnn. The magnetic moment of the cell is mainly contributed by the Mnj and 
Mnn atoms (/iMni = 2.21 /iB, /iMnn = 2.37 /ie), while the induced magnetic moments on 
Ge and C (/ice = —0.14 /xb, /Uc = —0.26 /ib) are much smaller and are antiparallel to 
the moments on Mn. The structure relaxations caused by interstitial atoms decrease the 
interatomic Mn-Mn distances (Tab. [I]), reducing the magnetic moments to fiuni = 1-99 /ib, 
fJ'Mnu = 2.14 /iB, /iGe = —0.13 /iB and fic = —0.21 /tB. However, the effect of hybridization 
on the decrease of the total magnetisation is much larger than the effect of relaxation. This 
conclusion follows from our electronic structure calculations for the relaxed MnsGesC system 
with the empty spheres at carbon positions (further referred to as MusGesVc). The calcu- 
lated magnetic moments on Mn agree well with the average saturated moment of 2.2 /iB/Mn, 
observed for the C-implanted MusGesCo.s films,- but deviate from those obtained for the 
sputtered samples, where the largest moment of 1.1 /iB/Mn was observed for MnsGesCo.rs^. 
This discrepancy could be ascribed to disorder effects and defects formation, which may 
differ in sputtered and implanted samples. 

In Tab. [T]we present some of the important calculated exchange constants for the relaxed 
and rigid MnsGes and MnsGesC systems (see Fig. [1] on the notation). Only the interac- 
tions between Mn spins are shown, as Mn-Ge interactions turned to be negligibly small. 
The first nearest neighbor Mnj-Mni interactions Ji (at the shortest distance of 2.526 A) are 
ferromagnetic (FM) and clearly dominate over the corresponding Mnj-Mnn and Mnn-Mun 
interactions, confirming the assumption about the different magnitude of exchange interac- 
tion for Mn atoms in different sublattices^. The dependence of the exchange parameters on 
the interatomic distances is stronger for the Mnn sublattice, while it is weaker for the Mnj 
one. As it follows from Tab. [H there is an antiferromagnetic (AFM) J3 interaction between 
the Mnn atoms within z = j and z = j planes (perpendicular to the c-axis), and FM J4 and 
Je interactions between Mnn belonging to different planes. The negative sign for J3 could 
be responsible for the small degree of non-coUinearity in MusGcs, observed in the previous 
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first-principles calculations^^. However, an increase of the corresponding Mnn-Mnn distance 
from 2.974 A to 3.057 A after structure relaxation leads to the change from AFM to FM 
J3 parameter, while the rest of the parameters are left almost unchanged. As a result the 
estimated value of the Tq = 400 K for the relaxed Mn^Ge-^ system is somewhat higher than 
for the rigid one of 320 K, the latter being in a very good agreement with the corresponding 
experimentally observed value of 304 K-. 

The results presented in Tab. |T] show that interstitial carbon, even without relaxation, 
significantly changes the value of J2, J3 and Jg, leaving the J4 and Ji interactions almost 
uneffected. In particular, the J3 interaction becomes strongly FM compared to the undoped 
MusGes. According to the DOS of MusGeaC in Fig. [21 this increase of the FM interac- 
tion can probably be explained by the strong p-d hybridization between Mnn and C states, 
which enhances the hopping of the 3d electrons from the partially occupied d-states of one 
Mnn atom to a 3d orbital of the neighboring Mnn. The Jg parameter, which in the doped 
compound corresponds to the 180° Mnn-C-Mnn interaction {dashed line on Fig. [1]), is es- 
sentially lower than the corresponding direct Mnn-Mnn exchange in the parent compound, 
pointing to the possible superexchange interaction between Mnn via carbon, as predicted by 
Goodenough-Anderson-Kanamori (GAK) rules.— According to these, 90° Mnn-C-Mun in- 
teractions should be FM, which is the case for J3 and J4. In this model interatomic exchange 
constants are strongly dependent on the orbital overlap, and hence sensitive to changes in 
the interatomic distances. It should be noted that GAK rules have been succesfuUy applied 
to explain the first-order FM-to-AFM transition in MnsGaCi-^ antiperovskite compounds,— 
in which the local environment of Mn atoms is similar to that in MusGesC. 

The transition temperature, calculated using the exchange parameters for the non-relaxed 
C-doped system, is essentially higher (Tc = 430 K) than in the undoped compound. This 
is a direct indication that the strongly enhanced FM stability is not related solely to the 
variation of the Mn-Mn interatomic distances. To distinguish the effect of relaxations from 
the electronic structure effect on the enhanced FM interaction in MusGesC, we compared 
the values of the exchange parameters in the MnsGesVc with those calculated for the parent 
compound. We can conclude that the essential decrease of the Mnn-Mnn distances due to 
carbon insertion reduces the FM exchange interaction, and has a particular influence on J3 
making it strongly AFM, thus destabilizing the coUinear magnetic structure in MnsGesVc 
(Tab. [I]). The relaxation makes J3 in MusGesC somewhat lower, but still high compared 
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to that in the pure compound. The J4 parameter is higher in relaxed MusGesC, pointing 
to the 90° FM superexchange, enhanced by the corresponding distance decrease. Although 
the strong hybridization between Mn 3d and C 2p orbitals leads to the contraction of the 
interatomic distances and lowers the longer-distance exchange constants, the FM Mn-Mn 
interactions prevail in the MusGesC, thus giving larger contribution to the Tq than the 
AFM ones. The calculated Tc=450 K is only 20 K larger than in the non-relaxed MnsGeaC, 
pointing to the small effect of the structural distortions on the Tq increase in MnsGes due to 
C-doping. The obtained value is in a very good agreement with the maximum value of the 
experimentally measured Tq of 445 Kr^. To summarize, the present studies clarified that 
the enhanced FM stability of C-doped MusGea compounds is attributed to the appearance 
of the 90° FM superexchange between Mn atoms mediated by carbon, while the structural 
distortions are found to be a secondary effect. 

I.S. acknowledges the CEA nanosciences program and P.M. the ESF SONS, Contract 
ERAS-CT-2003-980409, for partial funding The simulations were performed at the CEA 
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FIG. 1: (Color online) Crystal structure and exchange coupling scheme for the carbon-interstitial 
phase MnsGeaCx at x = 1 (full occupancy of carbon): projection on xy-plane (left), side view 
{right). Colored spheres denote Ge (grey), Mnj (dark blue), Mnn (light blue) and interstitial C 
(black) at the center of Mnn octahedra. Important exchange parameters for different Mn-Mn 
distances are shown in different colors. The black line indicates the edges of the unit cell. 
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FIG. 2: (Color online) Spin-resolved total and atom-projected DOS, calculated within the KKR 
method, for the MnsGea [thin line) and MnsGeaC (thick line) phases in FM configuration. Note 
the additional fine structure in the DOS of Mnn for MnsGeaC due to Mnn-C hybridization. 
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TABLE I: Calculated (within KKR) exchange constants Jij and Tq for the rigid and relaxed 
MnsGes and MnsGesC compounds. The corresponding Mn-Mn distances (in A) are shown 
in square brackets. Positive (negative) values characterize FM (AFM) coupling. Results for 
MnsGesVc are presented to distinguish between the structural and chemical effects of doping. 



Jij (mRy) 


MngGea (rigid) MngGes (relaxed) MnsGegC (rigid) MngGesC (relaxed) 


Mn^GeciVr 


T-Mni 


-Mni 


+2.14 [2.526] 


+2.09 [2.527] 


+2.24 [2.526] 


+2.34 [2.498] 


+2.24 [2.498] 


tMiii- 


-Mnii 


+0.59 [3.068] 


+0.59 [3.039] 


+0.90 [3.068] 


+0.69 [3.124] 


+0.40 [3.124] 


rMnii- 


-Mnii 


-0.15 [2.974] 


+0.13 [3.057] 


+0.84 [2.974] 


+0.51 [2.714] 


-1.80 [2.714] 


TMnir 
^4 


-Mnn 


+0.51 [3.055] 


+0.52 [3.082] 


+0.44 [3.055] 


+0.60 [2.948] 


+0.28 [2.948] 




-Mni 


-0.10 [4.148] 


-0.06 [4.148] 


-0.16 [4.148] 


-0.22 [4.118] 


-0.09 [4.118] 




-Mnn 


+0.69 [4.263] 


+0.65 [4.341] 


+0.07 [4.263] 


+0.22 [4.008] 


+0.76 [4.008] 




(K) 


320 


400 


430 


450 


non-collinear 



10 



